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INTRODUCTION 


This  report  summarizes  progress  during  the  first  twelve  months  of  an 
eighteen  month  study  concerning  imaging  from  transient  radar  signatures  as 
applied  to  cone-like  targets.  A  partial  data  set  has  been  acquired  for  a 
group  of  specific  shapes,  images  have  been  generated  using  the  data,  and 
the  characteristics  of  the  imaging  process  itself  for  two  separate 
approaches  have  been  investigated. 

This  research  is  an  outgrowth  of  past  studies  of  the  transient 
signatures  of  radar  targets  in  general,  and  of  methods  of  target  recognition 
using  transient  signatures  [1].  Past  efforts  had  already  established  several 
facts: 


Simple  relationships  exist  between  the  transient  radar  signature  and 
the  geometrical  properties.  For  imaging,  the  primary  relationship  used  is 
the  profile  function  (or  target  transverse  cross-sectional  area  vs  line  of 
sight  distance)  information  which  is  contained  in  the  ramp  response  signature. 
More  recent  imaging  research  is  also  making  use  of  target  resonance  properties, 
derived  from  the  decaying  portion  of  the  transient  response,  and  measured  as 
a  function  of  incident  E-field  polarization,  in  order  to  obtain  additional 
orientation  and  shape  information  [2],  Finally,  other  transient  response- 
derived  target  resonant  properties  are  proving  useful  in  studies  of  automated 
target  discrimination  [3J. 

The  primary  problem  in  applying  this  process  to  free-space  objects 
with  current  radar  technology  is  the  required  frequency  spectrum  of 
the  interrogating  radar  signal.  The  present  specification  is  for  a 
signal  covering  a  10:1  frequency  band,  whose  wavelengths  range  from 
approximately  four  times  to  0.4  times  the  target  maximum  dimension.  Thus, 
for  practical  targets  this  process  would  presently  require  construction 
of  an  ultra-broadband  HF  or  VHF  radar  system  to  acquire  the  signatures. 

This  is  the  principle  problem  to  be  attacked  in  the  proposed  research. 

Applying  these  same  techniques  to  arrive  at  target  identification 
processes  on  cone-like  targets  using  other  sources  of  target  scattering 
information  is  the  fundamental  goal  of  this  research. 

Two  broad  classes  of  target  scattering  information  are  being  con¬ 
sidered. 

1.  Normal  radar  systems,  with  interrogating 
frequencies  from  500  MHz  to  X  band _ 


In  any  radar  identification  system,  it  is  known  that  the  signature 
features  on  which  the  discrimination  process  is  based  are  crucial  to  the 
speed  and  error  performance  of  this  system.  It  is  postulated  that  transient 
response  analysis  and  target  imaging  analysis  which  makes  use  of  the  high 
frequencies  in  addition  to  the  low  frequency  data  used  previously  will  lead 
to  added  insight  into  future  selection  for  rapid,  accurate  discrimination 
using  high  frequencies  only. 
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2.  High-power  laser  radar  (LADAR) 

1  A  companion  research  effort  at  this  Laboratory  is  studying  non-linear 

effects  in  ordinary  materials  under  intense  laser  radiation.  This  effect 
could  be  used  to  "down-convert"  a  transient-modulated  laser  signal  to  produce 
the  equivalent  of  the  radar  transient  response  [4].  The  relationship  of 
this  laser-produced  transient  response  to  the  radar  impulse  response,  and 
consequent  imaging  and  target  identification  implications,  will  be  studied 
S  during  the  last  six  months  of  the  contract. 

Transient  response  information  is  fundamental  to  all  of  this  effort. 
Section  II  discusses  a  specific  set  of  target  shapes  which  were  chosen  for 
study,  and  the  experimental  and  analytical  work  involved  in  obtaining  the 
radar  scattering  data.  This  effort  is  concentrating  on  the  backscattered 
%  radar  cross-section  at  interrogating  directions  within  30°  of  nose-on 

incidence.  The  systems  being  used  for  data  measurement  are  described,  and 
analytic  techniques  used  in  the  calculations  are  briefly  discussed.  The 
characteristics  of  the  data  obtained  thus  far  are  discussed,  and  an  estimate 
of  the  accuracy  is  presented.  Some  of  the  data  are  also  compared  to 

other  measurements. 

I 

In  Section  III,  the  time-domain  ramp  responses  obtained  from  the 
data  are  discussed,  and  the  application  of  the  limiting-surface  imaging 
process  to  the  data  is  treated.  The  imaging  process  itself  is  currently 
being  examined  with  a  purpose  of  increased  automation.  The  current  state  of 
j  this  effort  is  also  presented. 

In  Section  IV,  our  research  on  inverse  scattering  concepts  is  discussed. 
The  Bojarski -Lewis  transform  is  reviewed.  It  is  shown  that  studies  of 
target  identification  with  limited  angular  and  frequency  information  are 
equivalent  to  analysis  of  windowing  using  the  Bojarski -Lewis  transform.  An 
extensive  set  of  results  are  presented  to  illustrate  the  effects  of 
!  windowing  on  some  simple  two  dimensional  shapes. 

Section  V  contains  a  summary  of  all  effort  on  this  contract  to  date, 
and  a  discussion  of  work  to  be  done  during  the  remainder  of  the  contract 
period.  Also  presented  are  some  possible  practical  applications  of  the 
I  research. 

Complete  frequency  domain  data  and  time  domain  waveforms  obtained  to 
date  are  included  in  Appendices  A  and  B. 
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II.  TRANSIENT  RESPONSE  OF  CONE-LIKE  OBJECTS 


A.  Prel iminary 

This  part  of  the  project  involved  the  collection  and  analysis  of 
data  to  be  used  in  the  imaging  of  the  six  cone-like  targets  shown  in 
Figure  1.  The  project  was  divided  into  two  main  tasks.  The  first 
task  was  the  collection  of  data  on  5  cm  long  models  of  the  six  targets 
using  the  Ten  Harmonic  Radar  system.  This  system  measures  both  the  phase 
and  amplitude  of  the  backscatter  from  the  set  of  targets  over  a  frequency 
ranging  from  1.1  GHz  to  11.0  GHz  in  ten  harmonically  related  steps.  With 
the  five  centimeter  long  cone-like  targets,  this  frequency  range  included 
data  from  the  Rayleigh  region  up  to  the  resonance  region  of  all  of  the 
targets. 

The  second  part  of  the  project  was  a  theoretical  study  of  the  radar 
backscattering  for  the  cone-like  targets.  The  results  of  the  theoretical 
study  allowed  the  experimental  data  to  be  critically  evaluated.  The  radar 
cross  section  has  been  determined  for  each  target  at  each  of  the  harmonically 
related  frequencies  and  for  each  aspect  angle  and  polarization.  In  most 
cases,  both  experimental  and  theoretical  results  are  available.  In  some 
cases,  however,  the  experimental  or  the  theoretical  results  are  unreliable 
and  only  one  of  the  two  methods  can  be  used. 

B.  Original  System  Checkout 

A  block  diagram  of  the  experimental  ten-harmonic  radar  backscatter 
system  [5,6]  is  shown  in  Figure  2.  In  this  system  a  single  L-band  microwave 
source  is  used  to  harmonically  generate  a  set  of  frequencies.  The  HP 
network  analyzer  receiver  local  oscillator  is  used  to  modulate  the  higher 
harmonics  (balanced  modulators),  and  also  to  modulate  (lower  sideband;  the 
L  band  signal  going  into  the  reference  port  of  the  HP  network  analyzer 
receiver  head.  The  result  is  that  the  receiver  head  of  the  HP  network 
analyzer  downconverts  each  of  the  harmonically  related  backscattered 
signals  to  the  same  IF  frequency.  The  output  of  the  HP  network  analyzer 
is  in  fact,  the  vector  sum  of  the  set  of  signals  being  received.  (In 
actual  practice,  the  transmitter  does  not  transmit  all  10  harmonics 
simultaneously.  Limitations  in  the  transmitter  output  multiplexer  require 
that  the  signals  be  transmitted  in  bands  of  4  or  5  frequencies  at  a  time.) 

A  photograph  of  the  bistatic  antennas  in  the  anechoic  chamber  with  an 
enlarged  (for  the  photograph)  sphere  target  on  the  foam  pedestal  is  shown 
in  Figure  3. 

The  processing  of  the  data  from  this  point  reduces  the  contribution  of 
the  background  and  antenna  coupling  in  the  received  signal  and  recovers 
the  individual  amplitude  and  phase  of  the  individual  harmonics.  The 
operational  procedure  is  to  move  the  target  slowly  toward  the  radar 
antennas  and  to  take  amplitude  and  phase  data  samples  under  computer  control. 
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Figure  1.  Outlines  of  targets  used  in  this  study 
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Figure  2.  Block  diagram  of  the  ten  harmonic 
radar  system. 
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Arrays  of  amplitude  and  phase  versus  position  data  for  the  targets  of 
interest  as  well  as  the  empty  support  pedestal  and  a  calibration  and  test 
sphere  are  stored  in  the  computer.  The  raw  data  are  Fourier  transformed 
to  yield  the  amplitude  and  phase  of  the  individual  harmonics.  Any  antenna 
coupling  or  scattering  from  stationary  objects  is  rejected  at  this  point 
as  a  DC  term.  The  target  measurements  then  have  the  empty  pedestal 
(background)  data  subtracted  and  the  final  result  is  calibrated  by  using 
the  calibration  sphere  data.  In  the  end,  we  have  the  amplitude  of  the 
targets  expressed  as  absolute  radar  cross  section  (in  cm.)  and  the  phase 
of  the  backscattered  return  relative  to  the  center  of  the  calibration  sphere. 
The  check  sphere  data  serves  as  an  estimate  of  the  overall  accuracy  of  the 
system. 

In  the  initial  stages  of  this  project,  a  number  of  preliminary  system 
checks  and  calibration  tests  were  performed.  It  was  necessary  to  confirm 
overall  system  linearity  so  that  the  Fourier  techniques  would  not  generate 
spurious  results.  It  was  also  necessary  to  confirm  the  accuracy  of  the 
target  positioning  system  so  that  the  subtractions  of  the  empty  pedestal 
(background)  data  would  be  valid.  The  sensitivity  and  signal  to  noise 
behavior  of  the  network  analyzer  also  needed  to  be  examined.  Any 
degradation  of  the  data  at  the  low  end  of  the  dynamic  range  of  the  target 
and  background  measurements  would  clearly  cause  severe  errors  in  the 
results.  The  ten  harmonic  range  was  finally  considered  to  be  functional, 
and  a  study  of  the  overall  system  accuracy  was  undertaken,  using  the  data 
on  the  check  sphere  target.  Some  of  the  results  are  shown  in  Figure  4. 

In  this  figure,  the  amplitude  and  phase  of  the  individual  cross-sections  for 
each  harmonically  related  frequency  are  plotted  using  a  polar  format.  As 
can  be  seen,  the  system  produces  data  with  accuracies  on  the  order  of  10%. 

(In  this  case,  accuracy  is  defined  as  the  magnitude  of  the  error  phasor 
divided  by  the  magnitude  of  the  known  theoretical  phasor.)  As  will  be 
shown  in  the  time  domain  section,  this  accuracy  is  sufficient  for  time 
domain  analysis  and  imaging  of  the  targets. 

C.  Data  Collection  Program 

Once  the  system  was  operating  with  confidence  the  collection  of  the 
experimental  data  was  begun.  As  was  mentioned  in  the  previous  section  of 
this  report,  the  actual  data  collection  was  carried  out  in  several  steps. 

The  sequence  of  data  collection  steps  is  outlined  below. 

(1)  Configure  radar  system  for  a  particular  harmonic  group, 

(i.e.,  1  ;2  through  5;6  through  10). 

(2)  Move  the  object  1/2  wavelength  (first  harmonic).  Position 
versus  amplitude  and  phase  (256  samples)  are  loaded  into  a  computer 
file.  ("Object"  in  this  case  may  be  the  empty  pedestal  (background), 
the  calibrate  sphere,  the  reference  sphere,  or  the  cone  target  of 
interest. ) 

(3)  Perform  step  2  for  the  background,  the  spheres,  and  the  cone 
target  at  aspect  angles  of  0,10,20,  and  30  degrees. 
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Figure  4.  Polar  diagram  of  error  distribution  of  the 
data  for  the  small  sphere  target. 
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(4)  Use  Fourier  transform  techniques  to  transform  the  position 
versus  amplitude  and  phase  computer  files  into  harmonic  number 
(frequency)  versus  amplitude  and  phase  computer  files. 

(5)  Subtract  the  background  data  from  the  sphere  and  target 
data  (phasor  calculation). 

(6)  Use  computed  values  of  amplitude  and  phase  from  the  reference 
sphere  for  the  computation  of  absolute  radar  cross  section  of  the 
check  sphere  and  the  cone  target. 

(7)  If  the  computed  values  of  the  cross  section  of  the  check 
sphere  are  within  tolerance  of  the  measured  values,  the  cone  target 
data  are  accepted  as  valid  and  the  absolute  radar  cross  section 

of  the  target  cone  (amplitude  and  phase)  is  stored  into  computer 
f  i  1 es . 

(8)  Once  all  targets  have  been  measured  in  this  way,  a  new 
harmonic  group  is  chosen  and  we  proceed  to  step  one  until  data 
for  all  harmonics  have  been  measured. 

The  final  result  of  these  measurements  is  a  set  of  computer  files  giving 
the  absolute  radar  cross  section  amplitude  and  phase  versus  harmonic 
number  (frequency)  for  each  of  the  targets  (including  the  check  spheres) 
at  each  aspect  angle.  These  results  are  tabulated  in  Appendex  A  of  this 
report.  Plots  of  part  of  these  data  are  given  in  Figures  5-10, 

The  data  shown  in  Figures  5-10  all  exhibit  certain  similar 
properties.  In  each  figure,  for  example,  the  transition  from  the  resonance 
region  to  the  Rayleigh  region  can  be  observed  as  a  transition  from  an 
oscillating  behavior  (as  a  function  of  frequency)  to  a  rapid  reduction 
of  amplitude  and  a  phase  which  approaches  180  degrees.  The  higher  frequency 
values  of  radar  cross  section  amplitude  of  1  to  3  cm  (1  to  1.7  cnr) 
corresponds  approximately  to  the  area  of  the  base  of  the  cone.  It  seems 
that  the  length  of  the  cone  is  not  an  important  determining  factor  of  the 
radar  cross  section  of  a  cone  target  when  the  aspect  angle  is  small.  We 
can  also  note  that  the  vertical  polarization  data  is  more  well  behaved 
than  the  horizontal  polarization  data.  Polarization  sense  is  defined  in 
this  case  with  reference  to  the  plane  containing  the  lines  from  the 
bistatic  antennas  to  the  target.  The  behavior  of  the  horizontal 
polarization  is  partly  caused  by  the  cone  targets  themselves  and  partly  caused 
by  the  increased  signal  to  noise  problems  encountered  in  the  horizontally 
polarized  bistatic  radar  antenna  configuration. 

A  final  remark  in  this  section  concerns  the  accuracy  of  the  overall 
data.  As  can  be  seen  from  the  figures,  the  spread  of  the  data  is  some¬ 
times  quite  large  for  the  various  aspect  angles,  or  from  harmonic  number 
to  harmonic  number.  These  large  variations  are  due  in  large  part  to  the 
system  errors.  It  is  possible  for  the  cone  targets  to  have  a  small  radar 
cross  section  at  the  same  time  that  the  check  targets  have  larger  cross 
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Figure  5a. 
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Amplitude  and  phase  data  on  cone  #1  as  a  function 
of  frequency,  (vertical  polarization) 
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Figure  5b.  Amplitude  and  phase  data  on  cone  #1  as  a  function 
of  frequency,  (horizontal  polarization) 
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Figure  6a. 
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Amplitude  and  phase  data  on  cone  #2  as  a  function 
of  frequency,  (vertical  polarization) 
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Figure  6b.  Amplitude  and  phase  data  on  cone  #2  as  a  function 
of  frequency,  (horizontal  polarization) 
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Figure  7a.  Amplitude  and  phase  data  on  cone  #3  as  a  function 
of  frequency,  (vertical  polarization) 
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Figure  7b.  Amplitude  and  nhase  data  on  cone  #3  as  a  function 
of  frequency,  (horizontal  polarization) 


AMPLITUDE  (CMj 

0.00  0.25  0.50  0.75  i.00 


X 

X 


01  23U56785  10 


HARMONIC  NUMBER 

Figure  8a.  Amplitude  and  phase  data  on  cone  #4  as  a  function 
of  frequency,  (vertical  polarization) 
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Figure  9a.  Amplitude  and  phase  data  on  cone  #5  as  a  function 
of  frequency,  (vertical  polarization) 
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Figure  9b. 
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sections.  This  can  cause  the  data  to  be  accepted  even  though  the  cone 
data  contain  large  errors.  Our  procedure  has  been  to  use  the  theoretical 
results  to  indicate  the  presence  of  errors  and  to  use  the  computed  values 
where  accurate  measurements  could  not  be  obtained.  This  procedure  will  be 
discussed  in  the  section  on  time  domain  computations. 

D.  Large  Target  Measurements 

In  order  to  obtain  data  in  the  high  frequency  regime  on  the  cone 
targets,  a  set  of  targets  eight  times  larger  than  the  original  set  was 
constructed.  It  was  not  possible  to  measure  the  cross  section  of  these 
targets  in  the  same  way  as  the  original  target  set,  however.  It  is 
expected  that  the  large  targets  are  much  more  sensitive  with  respect 
to  frequency  and  aspect  angle  than  the  small  targets.  Ten  measurements 
over  a  10  to  1  band  of  frequencies  will  not  be  sufficient  to  represent  the 
radar  behavior  of  the  large  targets. 

For  these  reasons,  it  was  decided  to  implement  a  swept  frequency 
radar  system.  In  such  a  system,  the  target  would  be  stationary,  and  the 
signal  to  noise  improvement  introduced  by  the  moving  target  data 
processing  system  would  not  be  available.  On  the  other  hand,  the  targets 
would  have  radar  cross  sections  approximately  64  times  larger  due  to  their 
8  times  larger  dimensions.  The  basic  concept  of  the  swept  frequency 
system  and  the  steps  in  the  measurement  of  the  data  are  shown  below. 

(1)  Measure  and  store  in  computer  files  amplitude  and  phase 
readings  as  the  frequency  is  stepped  over  a  band  of  frequency. 
Readings  are  taken  on  the  targets,  on  reference  and  check  spheres, 
and  on  the  background. 

(2)  Subtract  (phasor  subtraction)  the  background  data  from  the 
various  target  and  sphere  data  to  remove  the  errors  due  to  various 
antenna  coupling  and  unwanted  scatter  effects. 

(3)  Use  computed  values  of  the  reference  sphere  to  calibrate  the 
system  in  terms  of  absolute  radar  backscatter  cross  section 
(amplitude  and  phase). 

(4)  Use  data  on  the  check  sphere  to  confirm  the  overall  integrity 
of  the  system  and  the  data. 

The  major  problem  in  the  above  set  of  concepts  is  due  to  constraints  on 
the  frequency  repeatability  of  the  swept  frequency  system.  If  the  back¬ 
ground  is  to  be  subtracted  from  the  data  on  the  various  targets,  then  the 
frequency  for  each  step  in  the  frequency  sweep  must  be  repeatable  to  a 
high  degree.  The  constraint  is  that  the  amplitude  and  phase  of  the  back¬ 
ground  must  be  unchanged  when  the  target  is  introduced.  This  can  be 
assured  if  the  set  of  data  samples  is  taken  at  nearly  identical  frequencies 
for  tne  background  and  for  the  set  of  targets.  These  qualitative  comments 
were  studied  early  in  this  contract.  As  part  of  this  study,  measurements 
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were  made  of  the  relationship  between  phase  deviations  due  to  frequency 
deviations  and  analog  voltage  variations  due  to  frequency  deviations. 
Typical  results  are; 


Estimated  phase  accuracy  requirement 
Typical  resultant  frequency  resolution 
Resultant  analog  voltage  accuracy 
Absolute  voltage  accuracy 


=  -  5  deg. 

=  ±  0.003  GHz 
=  -  0.0088  volts/10  volts 
=  -  1  part  in  1100  (10  bits) 


The  analog  voltage  accuracy  discussed  in  the  above  list  is  the  voltage 
accuracy  that  would  be  needed  if  an  analog  voltage  was  used  to  control 
the  frequency  of  the  available  microwave  sweep  frequency  oscillator.  This 
required  voltage  accuracy  of  (typically)  1  part  in  1100  would  require  a 
design  accuracy  of  perhaps  1  part  in  2000  (11  bits).  This  stringent  analog 
voltage  stability  and  resolution  requirement  must  be  maintained  even 
though  the  relationship  between  voltage  and  frequency  is  dependent  on,  the 
sweep  oscillator  front  panel  dial  settings.  The  analog  system  would 
require  calibration  each  time  the  equipment  was  adjusted.  It  is  these 
limitations  which  resulted  in  the  decision  to  build  a  digital  interface 
between  the  digital  frequency  counter  and  the  computer.  Experiments 
indicated  that  the  drift  in  the  internal  oscillator  in  the  frequency 
counter  would  not  cause  significant  errors  in  the  system.  The  output  is 
inherently  of  high  resolution  (better  than  5  decimal  digits  (1  part  in 
10,000)  and  requires  no  calibration  at  all.  Implementation  of  the  digital 
system  has  been  completed  and  the  system  is  shown  in  block  diagram  form 
in  Figure  11.  A  computer  flow  diagram  of  the  software  that  was 
implemented  to  operate  the  digital  frequency  sweep  system  is  shown  in 
Figure  12. 


E.  Theoretical  Studies 


In  order  to  evaluate  the  measurements  in  an  orderly  fashion,  theoretical 
computations  of  the  radar  cross-sections  of  the  cone  targets  were  necessary. 
The  low-frequency  regime  is  represented  using  a  Rayleigh  calculation, 
while  the  high  frequency  regime  is  represented  using  the  Geometrical 
Theory  of  Diffraction  (GTD),  and  the  Equivalent  Current  Theory. 

The  Rayleigh  results  are  taken  from  Seigel  [7]  and  from  Crispin 
and  Siegel  [8],  The  Rayleigh  region  backscatter  cross  section  calculation 
derived  by  Seigel  is  given  below: 


a 


(i) 


where  f0  is  the  coefficient  of  the  leading  term  of  the  (poorly  known) 
series  expansion  in  powers  of  k0  of  the  exact  solution.  Seigel  gives  the 
final  solution  as 
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Figure  11.  Block  diagram  of  digital  frequency  sweep  radar  system 
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Figure  12.  Computer  flow  diagram  of  digital 
frequency  sweep  radar  system. 
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1  + 


4a  exp(^j-) 

tTFi 


where 


h  =  altitude  of  cone 
a  =  base  radius 


(2) 


(3) 


These  results  are  for  the  magnitude  of  the  radar  cross  section  only,  and 
no  phase  values  are  given.  The  phase  of  the  backscattered  signal  was 
approximated  in  this  study  by  taking  the  center  of  surface  area  (excluding 
the  base)  of  the  cone  targets  to  be  the  scattering  center.  (The  phase 
reference  zero  in  all  of  this  work  is  taken  at  the  base  of  the  cones.) 

There  are  a  large  number  of  references  on  the  Geometrical  Theory 
of  Diffraction  as  applied  to  the  finite  cone.  Most  of  the  references  are 
linked  to  the  paper  by  J.  8.  Keller  [9],  and  this  is  a  good  starting 
point  in  spite  of  the  number  of  corrections  to  this  paper  which  appear 
later  in  the  literature  (some  by  the  author  himself).  A  development  of 
the  Geometrical  Theory  of  Diffraction  (GTD)  is  beyond  the  scope  of  this 
report,  and  is  readily  available  in  the  literature.  Basically,  the  method 
generates  a  catalog  of  important  electromagnetic  ray  paths  and  considers 
the  behavior  of  the  primary  rays  of  the  radar  system.  The  diffraction  of 
rays  from  points  and  from  wedges  etc.  are  tabulated  and  combined  to  form  a 
representation  of  the  signal  backscattered  from  the  target.  In  the  studies 
presented  here,  the  cone  backscatter  was  represented  by  the  first  order  rays 
diffracted  directly  back  to  the  radar  by  the  edges  at  the  base  of  the  cone 
as  well  as  the  secondary  rays  diffracted  across  the  base  of  the  cone  and  then 
diffracted  again  back  to  the  observer.  The  equation  as  given  by 
J.  B.  Keller  [9]  is 


Eel (r,0)  +  Ee2(r,0)  = 

aftsin(^-)exp[ik(r+2acoty)] 

_ 


(cos  J  -  cos  jp) 


sin(l)exp(2ika-ij)  „  3^"2' 

- T7? -  *  (cos  rT  ~  cos  7i T' 

2  n  (»ka)1/2  n  2n  _ 


(4) 
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where  the  backscattering  cross  section 
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As  was  mentioned,  these  equations  consider  first  order  and  second  order 
backscatter. 

A  second  method  of  computing  the  radar  cross  section  of  the  cone 
targets  is  the  equivalent  current  theory  of  Burnside  and  Peters  [10], 

The  equivalent  current  referred  to  in  this  development  is  a  current  which 
radiates  the  same  fields  as  the  particular  target.  In  this  sense,  it  is 
not  a  target  surface  current.  The  diffraction  mechanism  considered  here 
is  similar  to  that  of  the  GTD  method  discussed  earlier  since  Burnside  and 
Peters  consider  the  same  first  order  and  second  order  diffraction  mechanisms. 
The  equations  which  describe  the  radar  cross  section  and  the  total  back- 
scattered  fields  are  given  below. 

2 

o=4Tra2  |-G(n,a)-  ^G(n,^-)’[Vg(n,2a,^)+Vg(n,2a,-^ +a)]  |  (6) 


with:  vR(n,if ,<>)  = 


expC^texpl-jkr) 


(27tkr) 


VT 


G(  n,  ip  ) 


and:  G(n,$  )  =  [~  sin(^)](cos^-  -  cos^-) 


*  -1 
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where:  <T  =  <ji  -  i|*0 

=  angle  from  cone  surface  to  observer 
ipQ  =  angle  from  cone  surface  to  source 
n  -  2(1  - 
a  =  cone  angle 
k  =  wave  number 
a  =  base  radius 
r  =  range  to  observer. 

Results  based  on  these  equations  will  be  given  shortly. 

Before  discussing  comparisons  between  these  computed  results  and 
the  experimental  data,  one  further  source  of  comparison  data  should  be 
discussed.  A  large  set  of  radar  cross  section  data  was  published  in 
1959  by  Keys  and  Primich  for  the  Canadian  Department  of  National  Defense 
[11].  The  cross-section  of  84  cones  was  measured  as  a  function  of  aspect 
angle  at  both  horizontal  and  vertical  polarization  (as  defined  earlier  in 
this  report).  The  cone  half  angles  ranged  from  4  to  20  degrees  and  the 
base  diameters  ranged  from  0.297  to  2.87  wavelengths  (in  14  steps).  Curves 
of  the  amplitude  of  the  radar  cross  section  as  a  function  of  the  aspect 
angle  are  available.  No  phase  data  is  available.  In  several  of  the 
comparisons  included  in  this  report,  values  from  this  Keys  and  Primich 
study  are  included.  It  should  be  remembered  that  the  accuracy  of  these 
values  is  on  the  order  of  i  2  dB. 

Comparisons  between  these  various  theories  and  measurements  and  the 
data  that  was  measured  in  this  program  will  be  discussed  next.  In  this 
part  of  the  report,  the  frequency  domain  results  and  comparisons  will  be 
given.  Time  domain  comparisons  will  be  given  later  in  this  report. 

Figures  13-16  show  the  comparisons  between  the  various  theoretical 
and  experimental  results.  In  these  figures,  the  experimental  cone  data 
are  all  taken  from  the  zero  degree  aspect  angle  data  sets.  In  fact,  the 
theoretical  results  shown  here  were  computed  for  monostatic  backscatter 
at  zero  degree  aspect  angle.  The  major  problem,  then,  is  to  reconcile  the 
monostatic  theoretical  values  with  the  measured  20  degree  bistatic  angles. 

It  can  be  shown  [8]  that  at  high  frequencies  the  monostatic  and  bi static 
results  are  equal.  Observation  of  our  comparison  between  theory  and  measure¬ 
ment  at  lower  frequencies  indicate  that  our  errors  are  small  at  low 
frequencies  as  well. 

The  general  behavior  of  the  set  of  comparisons  can  be  discussed  at 
this  point.  The  most  obvious  observation  is  that  there  is  often  a  region 
of  the  frequency  domain  data  where  none  of  these  theoretical  computations 
can  be  used  to  accurately  predict  the  behavior  of  the  cone  target.  In 
this  region  of  difficulty,  the  Rayleigh  computations  fail  because  the 
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Figure  15.  Comparison  between  theoretical  and  experimental 
results  for  the  amplitude  and  phase  versus 
frequency  for  target  number  3. 
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Figure  16.  Comparison  between  theoretical  and  experimental 
results  for  the  amplitude  and  phase  versus 
frequency  for  target  number  6. 
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frequency  is  too  high,  and  the  GTD  and  equivalent  current  theory  fail 
because  the  frequency  is  too  low.  It  is  in  this  regime  where  confidence  in 
the  measurements  is  very  important  since  there  is  no  independently  derived 
information  on  the  radar  cross  section.  Often,  this  region  where  the 
computations  fail  is  where  the  diameter  of  the  cone  is  from  1/50  to  1/4 
wavelengths.  In  general  it  can  be  seen  that  the  GTD  and  equivalent 
current  results  become  fairly  accurate  when  the  base  diameter  of  the  cone 
is  1/3  to  1/2  wavelength  or  greater.  This  implies  an  extension  of 
physical  optics  approximations  (of  which  the  GTD  and  equivalent  current  ~ 
solutions  are  special  types)  to  lower  frequencies  than  is  usually  the 
case.  « 

Conversion  of  these  frequency  domain  data  to  the  time  domain  will 
be  discussed  in  the  next  part  of  this  report. 

III.  RAMP  RESPONSES  AND  LIMITING  SURFACE  IMAGES 
A.  The  Ramp  Response  Waveform 

The  ramp  response  signature  is  the  basis  for  the  limiting  surface 
imaging  technique  used  in  this  study.  The  relationship  between  a  target's 
ramp  response  and  it's  geometrical  properties  was  suggested  by  Kenn-augh 
and  Moffatt  in  1965  [12].  The  basic  result  is  that 

FrU)  -  -  I  A(z)  (9) 


where  Fr(z)  is  the  time  domain  ramp  response  and  A(z)  is  the  target  cross 
sectional  area  versus  distance  along  the  direction  of  propagation  of"the 
incident  field.  This  relationship  was  first  exploited  in  a  radar  imaging 
technique  by  Young  in  1971  [13].  A  practical  radar  system  cannot  be 
used  to  measure  an  exact  ramp  response  due  to  the  bandwidth  requirements. 
The  alternative  used  is  to  measure  10  harmonically  related  samples  of  the 
targets  spectral  response.  These  samples  are  then  properly  weighted  and 
combined  to  construct  a  10-frequency  Fourier  series  approximation  of  the 
target's  ramp  response  signature.  The  expansion  for  the  periodic  ramp 
response  FR(t),  is  # 
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FR(t)  =  N  l 
K  n=l 


A 

— cos{nu  t  +  <j>  ) 
2  o  n' 

Fjt  n 


where  An  are  the  measured  amplitudes  of  the  spectral  response  samples,  $n 
are  the  measured  phases,  and  w0  =  (2tt ) (1 .08x10^)  rad/sec.  The  measurement 
system  is  documented  in  Section  II.  A  more  complete  discussion  of*the 
ramp  response  waveform  may  also  be  found  in  Reference  [13], 


B.  Ramp  Response  Processing 

The  ramp  response  derived  from  radar  cross  section  phase  and 
amplitude  versus  frequency  data  must  be  further  processed  to  yield 
useful  imaging  results.  The  modification  of  these  waveforms  is  based  upon 
known  geometrical  constraints.  The  ramp  response  was  shown  to  be  pro¬ 
portional  to  the  cross  sectional  area  versus  distance  of  the  target.  The 
response  should  therefore  be  zero  outside  the  main  region  of  the  response, 
i.e.,  the  region  corresponding  to  the  target.  However  the  response,  as 
constructed  above,  does  not  satisfy  this  because  it's  D.C.  term  was  not 
measured.  Previously,  the  D.C.  term  was  removed  by  shifting  the  waveform 
using  visual  estimation  and  the  CRT  display  and  light  pen.  However,  in  order 
to  study  error  associated  with  the  imaging  process  more  closely,  it  was 
necessary  to  systemize  this  procedure.  The  ramp  response  for  target  2,  a  4:1 
blunt  cone,  at  a  0°  look  angle  with  vertical  polarization  is  shown  in  Figure 
17a.  The  portion  of  the  waveform  ahead  of  the  main  response  is  averaged. 

The  entire  waveform  is  then  shifted  by  this  amount  as  in  Figure  17b. 

However  a  non-zero  response  still  exists  outside  of  the  target  region  due 
to  Gibbs  phenomena  ripple  introduced  by  bandlimiting.  These  portions  of 
the  waveform  are  set  to  zero  by  truncating  the  curve  at  the  points  where 
the  main  response  swings  negative.  This  is  shown  in  Figure  17c. 

When  imaging  with  multiple  look  angles  a  second  constraint  must 
also  be  satisfied.  The  volume  of  the  target  as  seen  from  each  look  angle 
must  be  consistent.  To  comply  with  this  requirement,  each  shifted  and 
truncated  waveform  for  a  particular  target  is  integrated  and  the  resulting 
volumes  are  averaged.  Each  waveform  is  then  scaled  so  that  the  volume  as 
seen  from  each  look  angle  is  equal  to  the  average  volume.  Figure  18  shows 
a  set  of  shifted  and  truncated  ramp  responses  for  target  2  for  0,10,20 
and  30°  look  angles  before  scaling.  The  volumes  for  each  are  indicated  and 
range  from  1.98  cm3  to  3.25  cm3.  The  average  is  2.58  cm3.  Figure  19  shows 
the  waveforms  after  scaling.  Once  a  set  of  waveforms  have  been  thus 
processed  they  are  ready  for  imaging. 

C.  Ramp  Responses  for  Cone-like  Targets 

Although  the  goal  of  radar  imaging  is  to  aid  in  target  discrimination, 
several  valuable  characteristics  can  be  obtained  directly  from  the  ramp 
response  signatures.  Figures  20a-f  are  scaled,  vertical  polarization, 
nose-on  ramp  responses  for  targets  1  through  6,  respectively.  The  actual 
target  profile  functions  are  included  for  comparison.  Due  to  a  significant 
phase  measurement  error  associated  with  the  first  harmonic,  it  was 
ncessary  to  alter  this  phase  value  before  constructing  several  of  the  ramp 
responses.  In  particular,  the  first  harmonic  phases  for  targets  1,3,4 
and  5  were  set  to  180°.  The  justification  for  this  stems  from  the  back- 
scatter  behavior  predicted  in  the  Rayleigh  region  as  discussed  in  Section 
II-E.  A  full  set  f  shifted,  unsealed  ramp  responses  for  each  target  are 
presented  in  Appendix  8.  Vertical  and  horizontal  polarization  ramp 
responses  are  included.  In  addition,  ramp  responses  are  included  with 
altered  and  unaltered  1st  harmonic  phases. 


Figure  17.  Ramp  response  modification. 
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Figure  20.  Vertical  polarization  ramp  responses  and  profile  functions 
i  from  nose  on  incidence  for  all  targets. 
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One  important  target  characteristic  which  can  easily  be  derived  from 
the  ramp  response  is  target  volume.  Table  III-l  contains  the  actual  volume 
for  each  target  along  with  the  volume  measured  from  the  ramp  response. 


TABLE  III-l 
TARGET  VOLUMES 


TARGET 

CALCULATED  VOLUME 
(cm3) 

MEASURED  VOLUME 
(cm3) 

VOLUME  ERROR 

1 

1.31 

1.56 

.16 

2 

2.04 

2.59 

.21 

3 

0.91 

0.98 

.07 

4 

1.41 

2.56 

.45 

5 

1.28 

2.20 

.42 

6 

3.64 

4.58 

.20 

In  each  case,  the  measured  volume  is  greater  than  the  actual  volume  due  to 
the  slow  decay  of  the  trailing  edge  of  the  waveforms.  However,  for  four 
of  the  targets  the  volume  error  is  less  than  .21.  These  four  targets, 
namely  targets  1,2,3,  and  6  could  be  discriminated  or  the  basis  of  volume 
alone.  The  disagreement  with  target  5  is  due  to  the  fact  that  this  target 
has  the  smallest  radar  cross  section  of  any  of  the  targets  and  is  too  small 
for  satisfactory  experimental  results. 

Most  of  the  targets  can  also  be  discriminated  from  the  slope  of  the 
ramp  response  leading  edge.  With  the  exception  of  target  1,  the  waveforms 
in  Figures  20a-f  can  be  ordered  according  to  slope  to  agree  with  the  slopes 
of  the  actual  profile  functions.  The  ramp  responses  in  order  of  increasing 
leading  edge  slope  are  5,  3,  4,  2,  and  6,  which  agrees  with  the  actual 
profile  functions.  Target  1  does  not  agree  because  when  the  waveform  was 
processed,  the  length  was  cut  off  along  the  0°  look  angle  more  than  along 
the  other  look  angles.  This  caused  the  0°  ramp  response  to  be  scaled  up 
excessively  during  the  scaling  process.  From  the  ramp  responses  for 
target  1  in  Appendix  B,  it  can  be  seen  that  the  slope  of  the  ramp  response 
leading  edge  corresponded  much  more  closely  to  the  actual  target  before 
scaling. 

Targets  1  and  4  are  identical  except  for  the  small  cylinder  on  the 
back  of  target  4.  As  can  be  seen  from  the  ramp  responses  in  Appendix  B, 
the  trailing  edges  of  the  target  4  ramp  responses  fall  off  more  gradually 
than  those  of  target  1.  The  cylinder  can  therefore  be  detected,  though 
not  in  detail . 

One  target  feature  is  not  detectable  in  the  ramp  responses.  Targets 
3,  5,  and  6  have  sharp  tips  and  back  edges,  while  targets  1,  2,  and  4  have 
blunt  tips  and  rounded  back  edges.  The  bandlimiting  of  the  ramp  response 
causes  rounding  off  of  all  sharp  discontinuities  on  the  targets,  thus 
overshadowing  these  details. 
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Figure  21  shows  the  unsealed,  horizontal  and  vertical  polarization 
ramp  responses  for  target  2.  As  the  look  angle  is  increased,  the 
horizontally  polarized  field  excites  a  resonant  mode  on  the  target.  The 
result  is  evident  in  the  ramp  responses.  This  resonant  mode  is  not 
predicted  by  physical  optics,  which  was  used  in  the  derivation  of  the  ramp 
response  approximation.  The  relationship  between  the  ramp  response  and 
target  profile  function  is  therefore  not  valid  when  the  incident  field  is 
polarized  along  the  major  axis  of  a  long,  narrow  target. 

0.  Description  of  Limiting  Surface  Imaging  Process 

The  limiting  surface  imaging  process  uses  one  or  more  profile 
functions  for  a  given  target  as  input.  The  idea  is  to  fit  a  surface  to 
the  set  of  profile  functions.  An  ellipsoidal  limiting  surface  is  used  for 
each  look  angle.  The  specification  of  a  series  of  input  parameters  is 
required.  Two  of  the  parameters,  e  and  <)>  specify  the  look  angle  for  a 
particular  profile  function.  The  remaining  parameters  are  in  general, 
unknown  and  require  a  priori  knowledge  of  the  target  geometry  or  iterative 
manipulation  to  achieve  the  best  image.  These  parameters  describe  the 
ellipsoidal  limiting  surfaces  to  be  used.  A  complete  set  of  these  parameters 
must  therefore  be  specified  for  each  look  angle.  A  straight  line  along 
which  the  centers  of  ellipses  are  to  be  aligned  must  be  specified.  This 
center  line  is  described  by  two  angles,  a  and  8.  a  is  the  angle  between  the 
center  line  and  the  look  angle  axis,  a  gives  the  rotation  of  the  center 
line  about  the  look  angle  axis  with  the  y  axis  as  the  8=0  reference. 

Gamma  describes  the  tilt  of  the  ellipses  about  the  center  line  and  is  also 
referenced  to  the  y  axis.  Each  of  these  angles  are  depicted  in  Figure  22. 

In  addition  the  aspect  ratio  for  the  ellipses  must  also  be  specified. 

Functionally  the  imaging  algorithm  can  be  described  as  follows.  For 
a  given  point  in  the  x-y  plane  the  program  iterates  through  values  of  z. 

For  each  z  value  the  area  of  an  ellipse  corresponding  to  each  look  angle 
is  calculated  which  passes  through  the  point  and  satisfies  the  parameters 
specified  for  that  look  angle.  If  each  of  the  areas  are  found  to  be 
less  than  the  values  of  the  corresponding  points  on  the  corresponding 
profile  functions,  then  the  point  is  assumed  to  lie  within  the  volume  of 
the  target.  The  points  with  minimum  and  maximum  z  values  which  are  found 
to  be  interior  to  the  target  are  chosen  to  lie  on  a  surface  contour  in 
the  y-z  plane.  A  complete  contour  is  swept  out  by  repeating  this  process 
while  iterating  over  y  and  holding  x  fixed.  Finally  a  complete  set  of 
contours  describing  the  target  surface  are  drawn  by  iterating  over  x. 

Portions  of  contours  which  should  be  visually  obscured  by  preceding 
contours  are  omitted  to  heighten  the  image's  3-dimensional  effect. 

The  correspondence  between  an  elliptical  cross  section  and  profile 
function  is  depicted  in  Figure  23.  The  look  angle  is  in  the  x-y  plane 
with  <t>=45°.  An  ellipse  perpendicular  to  the  look  angle  and  passing 
through  the  point  p  is  shown.  The  ellipse  has  area  a  which  is  also  the 
value  of  the  profile  function  at  the  corresponding  point  along  the  look 
angle  axis.  Therefore  the  point  p  lies  precisely  on  the  target's 
surface.  A  surface  contour  containing  p  is  also  indicated.  When  several 
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Figure  22.  Specification  of  angles  used  in  imaging  algorithm. 
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look  angles  are  used,  the  final  result  is  the  same  as  if  a  separate  image 
were  calculated  for  each  look  angle  and  their  intersection  taken  as  the 
actual  image. 


One  method  of  determining  image  quality  would  be  a  visual  comparison 
of  an  image  to  it's  corresponding  target.  However  a  more  objective  measure 
is  the  volume  error  between  the  input  and  output  profile  functions.  Volume 
error  is  defined  as: 


QO 

I  lAin<x)  -  A0Ut(x>l  dx 
.00 

eVOl.  ~ 

J  Ain(x)  dx 

—  oo 


(11) 


where  Ain(x)  is  the  ramp  response  or  input  profile  function  and  Aout(x) 
is  the  output  profile  function  for  the  image  along  a  given  look  angle. 

The  volume  error  is  calculated  for  each  look  angle  and  the  average  is 
used  as  a  basis  for  comparing  image  quality. 

Volume  error  is  strongly  influenced  by  the  specification  of  limiting 
surface  parameters  when  more  than  one  look  angle  is  used.  When  these 
parameters  are  not  correct  the  profile  functions  for  the  various  look 
angles  do  not  agree.  The  result  is  that  portions  of  the  target  are  cut 
off  in  the  image  and  the  entire  image  is  distorted.  However  when  only  one 
look  angle  is  used  there  are  no  other  profile  functions  to  agree  with, 
and  the  volume  error  becomes  virtually  independent  of  the  limiting  surface 
parameters.  Since  the  algorithm  essentially  calculates  elliptical 
cross  sections  whose  areas  correspond  to  specific  points  on  the  input 
profile  functions  it  would  seem  that  there  should  be  perfect  agreement 
between  input  and  output  profile  functions  in  this  case.  However  a 
secondary  source  of  volume  error  becomes  apparent.  Due  to  the  discrete 
nature  of  the  image  surface  and  hence  the  image  cross  sections,  quantization 
error  is  introduced  into  the  output  profile  function.  Strictly  speaking, 
quantization  error  should  have  zero  rms  value.  However  because  of  the  way 
in  which  the  limiting  surface  algorithm  operates,  the  output  profile  function 
is  always  less  than  the  input  ramp  response.  The  quantization  error  therefore 
has  a  negative  bias.  Figure  24  depicts  this  quantization  error.  The  input 
and  output  profile  functions  for  an  image  of  target  2  using  a  0  degree  look 
angle  are  shown.  Quantization  error  in  the  output  profile  is  evident.  The 
volume  error  is  0.11.  This  indicates  that  for  several  look  angles,  even  if 
all  parameters  were  known  and  the  target  had  elliptical  cross  sections 
throughout  it's  entire  volume  from  all  look  angles,  such  as  an  ellipsoid, 
volume  errors  less  than  the  order  of  0.1  could  not  be  expected. 
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For  general  targets  that  are  not  comprised  of  a  few  simple  shapes, 
the  calculation  of  the  limiting  surface  parameters  discussed  earlier  is 
not  possible.  At  best  they  can  be  only  roughly  approximated.  However  for 
the  class  of  targets  in  this  study,  i.e.,  conical  shapes,  the  parameters 
can  be  closely  approximated.  This  is  because  the  cone  cross  sections  are 
ellipses  along  a  limited  range  of  look  angles.  If  ip  is  the  generating 
angle  of  an  infinite  cone  and  <j>  is  the  look  angle  from  the  tip  of  the 
cone  with  respect  to  the  cone's  axis  of  symmetry,  then  the  cross  section 
perpendicular  to  the  look  angle  will  be  an  ellipse  providing  0< |(j> | <90°-xp. 
For  e=0  the  cross  section  will  be  a  circle.  This  look  angle  constraint 
is  satisfied  in  each  of  our  cases.  However  since  our  cones  are  finite, 
the  elliptical  cross  sections  at  the  back  of  the  cones  are  partially 
truncated.  The  calculated  parameters  will  agree  precisely  at  the  tip  of 
the  cone  but  will  cause  a  certain  amount  of  distortion  near  the  base. 

The  eccentricity,  e,  of  an  ellipse  indicates  it's  degree  of 
departure  from  circularity.  Eccentricity  varies  between  0  for  a  circle 
and  1  for  a  straight  line.  For  an  elliptical  conic  section,  eccentricity 
can  be  expressed  as  [14], 

„  _  cos(90-») 


However  eccentricity  can  also  be  expressed  in  terms  of  the  ellipse's  major 
axis  half  length  a,  and  the  distance,  c,  from  the  ellipse's  center  to  one 
of  the  foci , 
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A  third  relation  relating  a,  c,  and  the  ellipse's  minor  axis  half  length,  b 
is, 

c2  =  a2  -  b2  .  (14) 

Combining  Equations  (12),  (13),  and  (14)  above  we  have. 


where  ^  is  the  aspect  ratio  for  the  ellipse.  For  target  2,  tp  is  6°.  The 
aspect  ratios  for  the  various  look  angles  are  shown  in  Table  1 1 1- 2 


TABLE  1 1 1-2 

ELLIPTICAL  CROSS-SECTION  ASPECT  RATIOS  VS. 
LOOK  ANGLE  FOR  TARGET  2 


LOOK  ANGLE 

ASPECT  RATIO 

4> 

a/b 

0° 

1.0 

10° 

1.06 

30° 

1.16 

The  line  containing  the  center  of  the  elliptical  cross  sections  is 
shown  in  Figure  25.  For  the  range  of  <j>'s  and  <|/s  that  we  are  interested 


Figure  25.  Location  of  elliptical  cross  sections  for  a  cone. 
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in,  the  angle  <5  between  the  center  line  and  the  cone's  axis  of  symmetry 
is  very  small.  Therefore  a  is  approximately  the  same  as  the  look  angle. 
Due  to  the  orientation  of  the  target  along  the  x  axis  and  the  fact  that 
the  look  angles  all  lie  in  the  x-y  plane,  6=180°  and  y=0°  in  all  cases. 


Figure  26  is  an  image  of  target  2  using  the  0  and  30  degree  look 


Figure  26.  Image  of  Target  2  generated  from  ramp  responses  from 
0°  and  30°  look  angles  with  optimum  parameters. 


angle  profile  functions  from  Figure  19.  The  undulation  in  the  surface  is 
due  to  the  ripple  in  the  profile  functions  resulting  from  the  10  harmonic 
ramp  response  approximation.  Distortion  can  be  seen  near  the  base  of  the 
image  due  to  the  truncated  elliptical  cross  sections  of  the  cone  along 
the  30°  look  angle.  Figure  27  shows  the  comparison  of  the  input  and 


Figure  27.  Input  and  output  profile  functions  for  image  in  Figure  15. 
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output  profile  functions  for  the  image.  The  back  end  of  the  output  profile 
function  along  the  zero  degree  look  angle  is  seen  to  be  significantly  cut 
off  in  comparison  with  the  input  profile  function 

Some  images  generated  for  the  other  targets  are  shown  in  Figures 
28  through  32.  Each  of  these  images  was  generated  using  the  0°  look  angle 
(nose-on)  ramp  response  waveform.  The  volume  error  estimate  is  shown  along 
with  the  image.  In  general,  all  images  have  those  features  just  discussed 
in  the  case  of  target  2. 

F.  Imaging  Dependence  on  Limiting  Surface  Parameters 

As  was  pointed  out  earlier,  the  limiting  surface  parameters  for  a 
given  target  are  initially  unknown.  However  to  obtain  an  accurate  image 
they  must  be  correctly  determined.  One  possible  approach  would  be  to  try 
various  combinations  based  upon  reasonable  guesses  until  a  suitable  image 
resulted.  This  approach  is  not  satisfactory  for  unknown  targets.  An 
alternative  is  to  study  the  agreement  between  an  image  and  it's  input  profile 
functions.  The  volume  error  is  a  convenient  measure  of  this  agreement.  The 
volume  error  can  be  observed  as  a  function  of  one  or  more  limiting  surface 
parameters.  The  goal  will  be  to  obtain  quantitative  information  on  the 
relation  between  volume  error  and  the  limiting  surface  parameters. 

To  investigate  this  relationship,  images  were  generated  using  the 
0  and  30  degree  look  angle  profile  functions  for  target  2.  These  are 
shown  in  Figure  19.  In  each  case  correct  values  of  parameters  for 
the  0  degree  look  angle  limiting  surface  were  used.  Parameters  for  the 
30°  look  angle  limiting  surface  were  iterated  one  at  a  time  while  setting 
the  others  at  their  correct  values.  For  the  0  degree  look  angle  the 
parameters  are;  a=0,  6=0,  y=0  and  aspect  ratio  =  1.0.  For  the  30  degree 
look  angle  a=30°,  6=180°,  y=0°  and  aspect  ratio  =  1.16.  Figure  33 
shows  the  variation  in  volume  error  as  a  varies  from  27  to  36  degrees. 

The  minimum  error  is  0.215  and  occurs  at  a=30  degrees.  Figure  34 
demonstrates  how  volume  error  is  introduced  by  incorrect  values  of  a. 

Outlines  of  images  generated  independently  from  9  and  30°  profile  functions 
are  shown.  With  a=30°  the  images  are  properly  aligned.  A  second  case 
is  shown  with  a=45°.  Since  the  composite  image  is  the  intersection  of  the 
two  independent  images  it  is  seen  that  the  front  of  the  cone  will  be  clipped 
off  in  the  image.  Figure  35  is  the  resulting  image  and  may  be  compared 
with  Figure  26.  The  volume  error  is  0.55.  the  profile  functions  for 
the  misaligned  case  are  shown  in  Figure  36.  Figure  37  shows  the 
variation  in  volume  error  with  aspect  ratio.  The  aspect  ratio  is  incre¬ 
mented  from  0.25  to  3  with  a  minimum  volume  error  of  0.218  occurring  at 
1.125  which  is  very  close  to  the  calculated  value  of  1.16.  Figure  38 
is  a  plot  of  volume  error  versus  y.  The  curve  is  symmetric  about  y=0°  so 
only  the  range  from  0°  to  90°  is  shown.  The  volume  error  is  a  minimum  of 
0.215  at  y=0°.  The  volume  error  tends  to  decrease  as  y  approaches  zero 
although  the  curve  does  not  decrease  monotonical ly.  It  should  be  noted 
however  that  the  fluctuations  or  ripples  in  the  curve  are  comparatively 
small,  less  than  0.01  and  the  entire  range  of  volume  error  is  much  less 
than  that  seen  with  the  other  parameters. 
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VOLUME  ERROR  =  0.340 


Figure  28.  Image  of  target  1,  generated  from  ramp  responses 
at  0°  and  30°  look  angles. 


VOLUME  ERROR  =  0.444 


Figure  29.  Image  of  target  3,  generated  from  ramp  responses 
at  0°  and  30°  look  angles. 


Figure  30.  Image  of  target  4,  generated  from  ramp  responses 
at  0°  and  30°  look  angles. 


VOLUME  ERROR  =  0.326 


Figure  31.  Image  of  target  5,  generated  from  ramp  responses 
0°  and  30°  look  angles. 


Figure  32.  Image  of  target  6,  generated  from  ramp  responses 
0°  and  30°  look  angles. 


Figure  35.  Image  of  target  2  generated  from  ramp  responses  from  0° 
and  30°  look  angles  with  alpha=45°  along  30°  look  angle. 
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Figure  36.  Input  and  output  profile  functions  for  image  in  Figure  35. 
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Figure  38.  Volume  error  versus  y. 


With  each  of  the  three  parameters  studied  the  volume  error  was  minimized 
as  the  correct  parameter  value  was  reached.  In  two  of  the  cases  the  volume 
error  decreased  monotonically  towards  the  minimum  point  while  in  the 
third  case  there  was  a  trend  towards  decreasing  error  although  some  ripple 
was  present.  This  suggests  that  an  automated  iterative  technique  which 
seeks  to  minimize  volume  error  using  a  standard  approach  like  Newton's 
method  may  be  appropriate  for  determining  the  unknown  parameters.  However, 
volume  error  dependence  on  the  limiting  surface  parameters  must  be  studied 
in  greater  depth  before  this  is  attempted.  In  particular  volume  error 
must  be  observed  as  several  parameters  stray  from  their  optimum  values. 

This  will  be  the  focus  of  the  limiting  surface  imaging  study  in  the  near 
future. 

G.  A  Possible  Approach  for  the  Study  of 
Image  Confidence 

Image  quality  has  been  studied  by  comparing  input  profile  functions 
to  the  image  profile  functions.  However  this  may  not  be  a  complete  measure 
of  the  correspondence  between  a  target  and  it's  resulting  image.  The 
question  that  remains  is  whether  or  not  a  set  of  profile  functions  uniquely 
define  a  target.  To  study  this  problem  a  tv/o  dimensional  example  will  be 


For  one  look  angle  an  infinite  number  of  objects  can  be  found  to 
satisfy  a  given  profile  function.  This  can  be  demonstrated  by  slicing  an 
object  perpendicular  to  a  look  angle.  The  object  may  be  modified  by  dis¬ 
placing  the  slices  perpendicular  to  the  look  angle  in  any  arbitrary  fashion. 
An  example  is  shown  in  Figure  39.  The  profile  function  along  the  look  angle 
remains  unchanged. 

The  problem  is  constrained  when  a  second  look  angle  is  added.  It 
must  be  determined  if  it  is  possible  to  displace  the  slices  along  one  look 
angle  in  such  a  way  as  to  maintain  the  original  profile  function  along  a 
second  look  angle  as  well.  The  object  considered  will  be  a  rectangle. 

The  radar  resolution  will  be  assumed  to  be  t.  The  slices  will  therefore 
have  thickness  t  and  only  displacements  of  multiples  of  t  will  be  allowed. 

A  rectangle  sliced  along  an  off  axis  look  angle  is  shown  in  Figure  40. 

For  t>w  sin<i>,  where  w  is  the  rectangle  width  and  <t>  is  the  look  angle  from 
the  x  axis,  each  slice  will  have  3  distinct  regions  as  seen  along  the 
x  axis.  Assuming  the  back  tip  of  the  slice  to  be  at  x=0,  the  cross 
sectional  area  versus  x  for  the  slice  is: 

0  <  x  <  t/cos<f> 
t/cos 4>  <  x  <  wtan<p 

(16) 

wtan<j>  <  x  <  wtan<p  +  - 

COS<J> 

otherwi se 


A(x)  = 


x/tan(p 

t/sin<|> 

_t _  _  (x-wtan<j>) 

sin<p  ”  tan<f> 

0 


A(x)  is  plotted  in  Figure  41. 


In  order  to  avoid  having  slices  scattered  throughout  all  of  space  we  will 
require  that  some  part  of  each  slice  at  least  touch  the  original  volume 
of  the  rectangle.  The  maximum  displacement  of  slices  perpendicular  to 
the  look  angle  axis  determines  how  many  slices  ahead  one  must  look  for 
area  contribution  at  a  given  point  along  the  x  axis.  For  example  in 
Figure  42  it  is  seen  that  point  1  on  the  x  axis  may  receive  an  area 
contribution  from  any  slice  along  the  look  angle  axis  up  to  slice  6.  The 
area  contribution  at  a  given  point  on  the  x  axis  due  to  a  slice  along  the 
look  angle  axis  depends  on  the  actual  amount  by  which  the  slice  has  been 
displaced.  Assume  slice  j  along  the  look  angle  axis  is  displaced  by  nt. 

The  distance  from  the  tip  of  the  slice  to  point  i  along  the  x  axis  is 
jtcos<p-nt  sin<f>-t(i-l/2).  This  distance  is  then  inserted  into  the  formula 
for  A(x)  above  as  x.  A(x)  for  each  contributing  slice  must  be  found  for 
a  given  point  on  the  x  axis.  At  each  point,  one  must  look  ahead  by  a 
number  of  slices  equal  to  the  next  highest  integer  value  of  2wtan<p/t  cos<t>=N. 

i+N-l 

The  general  form  of  the  area  at  each  point  will  be  a.  =  \ 

1  J=i 


A..(x)  where 


Figure  39.  Profile  function  for  a  rectangular  solid 
before  and  after  mod  if i cation. 


Figure  40.  Rectangle  sliced  perpendicular  to  look  angle 
with  one  slice  displaced. 


Figure  41.  Cross  sectional  area  versus  distance  for  a 
slice  with  unit  thickness. 


Figure  42.  Area  contribution  at  point  1  along  7  direction 
due  to  slices  1-6  along  look  angle  direction. 


A-jj(x)  is  the  area  contribution  of  slice  j  at  point  i  along  the  x  axis. 

If  the  rectangle  has  m  slices  along  the  x  axis  then  there  will  be  m  such 
equations.  The  problem  is  to  solve  these  equations  for  the  combinations 
of  displacements  which  yield  the  same  a-j  as  the  original  rectangle. 
Unfortunately,  due  to  the  form  of  A(x)  for  each  slice,  the  relationship 
between  the  a^  and  the  slice  displacements  is  non-linear  and  a  closed  form 
solution  is  not  apparent.  However  one  observation  can  be  made.  As  the 
number  of  look  angles  increase  the  problem  becomes  more  constrained  and 
the  probability  of  several  objects  satisfying  the  same  set  of  profile 
functions  decreases. 


IV.  HIGH  FREQUENCY  IMAGING  USING  THE 

BOJARSKI  IDENTITY 

A.  Introduction 

As  we  indicated  in  the  proposal  for  our  present  research,  we  would 
endeavor  to  study  techniques  to  produce  target  images  based  on  the  backscatter 
data.  Numerous  researchers  have  reported  in  this  field  for  years  with 
several  approaches  to  the  same  basic  problem;  that  of  relating  the  back¬ 
scatter  to  the  phenomena  that  gives  rise  to  the  fields  and  thereby 
determine  a  satisfactory  image.  As  is  often  the  case  each  approach  is 
based  on  fundamental  assumptions  about  the  shape  of  the  body,  material 
properties  of  the  body  and  the  frequency  range  of  the  interrogating  signal. 

In  our  contribution  we  have  utilized  the  Bojarski  identity  [15]  which 
is  valid  in  the  physical  optics  limit,  i.e.,  that  the  currents  induced 
on  the  body  are  2nxTT. 

In  Section  IV-B  we  will  introduce  the  Bojarski  identity. 

In  Sections  IV-C  through  IV-E  we  will  generate  data  for  trial  of  the 
identity.  First  a  one  dimensional  characteristic  function  with  the  effects 
of  a  one  dimensional  bandpass  filter  to  simulate  the  frequency  constraint 
of  a  radar  system  at  one  aspect,  then  two  dimensional  characteristic 
functions  are  generated  for  triangular,  square  and  circular  shapes  and 
the  corresponding  spectrum  is  windowed  in  aspect  and  in  frequency  to  see 
the  resulting  distortion  in  the  characteristic  function  space.  These 
results  will  not  only  give  us  a  "feel"  for  the  images  to  expect  with 
finite  backscatter  data  but  they  also  demonstrate  the  Gibbs  distortion 
due  to  the  window  itself  and  through  the  proper  use  of  spectral  tapering 
these  effects  can  be  minimized. 

Next,  in  Section  IV-F  a  study  is  made  of  the  possibility  of  extending 
the  frequency-aspect  window  in  order  to  improve  the  images.  Specific 
attention  is  given  to  observing  the  nature  of  the  angular  variations  of  the 
scatter  to  analytically  extend  the  window  or  to  determine  which  features  are 
of  greater  significance. 


Finally  suggestions  are  made  for  further  study  and  insights  in  image 
improvement  are  discussed. 

B.  The  Bojarski  Identity 

In  a  series  of  reports  Norbert  Bojarski  [15,16,17,18,19]  has  formulated 
a  relationship  between  the  monostatic  scattered  far  field  cross  section 
of  perfect  conductors  and  the  geometry  of  the  conductors.  The  results  are 
based  on  the  physical  optics  approximation. 

If  the  scatterer  is  expressed  in  three  dimensions  by  its  characteristic 
function  y(3T)  where. 


) 1  inside  the  body 
Y'x'  \0  outside  the  body 


(17) 


L 

and  if  the  scattered  far  field  cross  section  p(p)  is  used  to  define  F(p] 
where. 


r(p)  =  /CT  £L J) 

I  \c 


(18) 


with  p(fT)  =  the  field  cross  section 

V  =  wavenumber  propagation  vector  in  the  direction 
of  aspect; 

then  Bojarski  has  shown  that 


'(x)  - 


(2tt)' 


oo 

J  r(p)  e1^"* 


x  ,3— 
d  x 


(19) 


which  is  just  a  3-dimensional  inverse  Fourier  transform  of  r(p). 

C.  Trial  Results 

At  the  outset  this  identity  suggests  that  we  need  £jp)  to  be  measured 
at  all  frequencies  and  aspects  to  exactly  reconstruct  v(x)  and  thus  the 
target.  But  there  are  many  factors  that  do  not  necessitate  knowledge  of 
the  complete  backscatter  spectrum  [20,21,22].  Tin's  is  seen  to  be  true  in 
mathematically  related  fields  where,  e.g.,  intelligible  communication  signals 
can  be  heavily  bandlimited  and  yet  the  information  can  be  satisfactorily 
recovered  within  some  tolerable  level. 


D.  One  Dimension 

Whenever  a  function's  spectral  components  are  bandlimited,  often 
termed  windowing,  then  the  inverse  transform  of  the  windowed  spectrum  will 
produce  a  likeness  of  the  original  function.  Certain  features  will  be 
preserved  depending  on  what  portion  of  the  infinite  spectrum  is  windowed. 

In  the  practical  radar  sense  only  a  finite  frequency  band  is  used  thus  it 
is  instructive  to  see  the  distortion  introduced  in  the  corresponding 
characteristic  function.  Figure  43  shows  a  one  aspect  characteristic 
function  (where  y(x)  represents  the  linear  extent  of  the  target  along  that 
aspect),  its  unwindowed  and  windowed  spectrum,  and  the  resulting  image  after 
windowing.  It  can  be  seen  that  impulses  occur  at  the  discontinuities.  Thus 
the  scatterers  front  and  back  edge  could  still  be  determined.  The  advantages 
of  the  preserved  edges  are  twofold;  first,  we  must  rely  on  this  result  in 
bandlimited  radar  systems  and  secondly,  the  low  frequency  components  in 
real  systems  do  not  obey  the  physical  optics  approximation  so  we  could  not 
use  them  with  extreme  accuracy  anyway.  And  if  we  can  determine  edges  in 
3-dimensions,  i.e.,  the  surface  of  the  scatterer,  then  we  have  exactly 
specified  the  shape. 

E.  Two  Dimensions 

We  have  elected  to  study  the  application  of  Bojarski's  identity  in 
two  dimensions  because  of  the  relative  ease  of  display  and  we  can  infer 
results  to  the  three  dimensional  case. 

The  method  of  display  of  the  two  dimensional  function  is  a  gray  level 
plot.  Each  array  of  values  in  either  x  or  p  space  is  quantized  into 
14  levels  and  each  level  is  assigned  a  gray  density  in  the  form  of  over¬ 
striking  characters  on  the  line  printer.  (A  blank  for  level  1,  and 
the  characters  U,T,X,$  overstruck  for  level  14.)  Additionally  the  array 
is  interpolated  in  the  V  dimension  to  expand  the  display  so  that  NxN  point 
arrays  appear  square  visually. 

Evaluation  of  triangular,  square  and 

circular  characteristic  function's 

In  this  section  we  demonstrate  the  results  of  applying  the  identity 
to  three  shapes.  An  equilateral  triangle  is  shown  in  Figure  44a.  The 
y  axis  dimension  is  11  spatial  samples  high  and  the  x  axis  dimension  is 
21  spatial  samples  wide.  The  actual  dimensions  are  omitted  because  the 
examples  demonstrate  principles  and  are  not  intended  for  a  specific 
application.  (Recall  that  the  y  axis  values  are  interpolated  by  the  gray 
level  plotting  routine  thus  more  than  11  samples  appear.)  The  “black" 
areas  represent  a  magnitude  of  1,  otherwise,  the  magnitude  is  0.  Thus  we 
see  a  "top"  view  of  the  characteristic  function  in  a  plane.  The  square 
characteristic  function  appears  in  Figure  44b  and  both  dimensions  are  11 
spatial  samples  long.  Finally  a  circular  characteristic  function  appears 
in  Figure  44c  with  a  diameter  of  11  spatial  samples.  After  application  of 
a  two-dimensional  FFT  the  magnitude  of  the  spectra  for  all  three  shapes 
appears  in  Figures  45,  46  and  47. 
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Figure  44c.  Characteristic  function  of  a  circular 
scatterer  in  two  dimensions. 


Two  dimensional  scattering  from 
the  triangular  target. 


It  is  noted  that  the  different  geometries  have  characteristic  back- 
scatter  variations.  The  square  and  triangle  have  large  return  "lobes" 
which  are  perpendicular  to  the  straight  edges.  These  correspond  to  the 
"speculars"  in  normal  radar  terminology.  The  resonances  of  the  circle  are 
reflected  in  oscillations  which  are  radially  directed  and  are  symmetric  in 
magnitude  about  the  origin.  Thus,  the  effects  of  degradation  due  to 
windowing  will  depend  upon  both  the  geometry  of  the  scatterer  and  the  spectral 
angular  location  of  the  windows.  This  may  be  reflected  in  two  ways.  One; 
by  the  law  of  conservation  of  energy,  if  the  geometry  gives  rise  to  lobes 
in  the  spectrum  then,  for  a  given  window  size,  varying  aspect  locations  will 
produce  varying  amounts  of  energy.  Thus  if  the  window  is  located  in  the 
vicinity  of  a  null  less  energy  is  available  for  image  reconstruction  than 
near  a  maximum.  Secondly,  lobes  are  perpendicular  to  surfaces  along  major 
dimensions  and  nulls  are  located  at  aspects  where  there  are  few  or  no 
corresponding  perpendicular  surfaces.  Therefore  the  most  significant 
measurements  would  be  taken  when  the  aspect  angle  of  the  receiver  lies  in 
the  region  of  major  lobes.  This  is  evidenced  by  the  test  results  in  the 
next  three  sections. 

Additionally,  we  know  that  the  characteristic  function  in  x”  space 
is  a  constant  (=1)  within  the  target  and  has  a  jump  discontinuity  at  the 
boundary.  Thus,  the  gross  target  volume  information  is  primarily  concentrated 
in  the  low  frequency  backscatter  and  the  location  of  surfaces  (the  profile) 
is  concentrated  in  the  higher  frequencies.  Therefore,  if  we  retain  only 
higher  frequency  information  we  may  still  reconstruct  the  surface  of  the 
scatterer  for  identification.  This  is  contrary  to  the  results  of  Perry's 
analysis  [23]. 

In  the  following  sections  we  consider  the  effects  of  finite  apertures 
or  windows  on  the  reconstructed  images. 

Frequency  windowing 

Since  practical  radar  systems  operate  in  a  band  of  frequencies  we 
have  simulated  this  by  applying  a  frequency  band  window  at  all  aspects. 

The  basic  window  appears  in  Figure  48  and  has  magnitude  1  within  the 
upper  and  lower  cutoff  and  0  in  the  stop  band,  this  window  is  then 
multiplied  with  the  scattering  spectrum  to  simulate  the  measurements  of 
a  bandlimited  radar.  In  Case  I  we  applied  a  window  with  a  lower  cutoff 
at  the  10-th  sample  point  and  the  upper  cutoff  at  the  30-th  sample  point. 

(Our  spectrum  is  32  points  long  in  the  positive  Px  and  positive  Py 
directions.)  The  resulting  images,  for  the  triangular,  square  and  circular 
shapes  appear  in  Figures  49  to  51.  Since  a  blank  represents  a  minimum  in 
the  gray  scale  plot  we  see  that  the  minimum  is  the  under  shoot,  due  to  Gibbs 
phenomenon,  at  the  edge  of  the  targets.  A  maximum  overshoot  also  occurs 
at  the  edges.  Therefore  the  boundary  information  is  still  preserved.  We 
applied  a  threshold  filter  at  a  level  approximately  half  way  between  the 
minimum  and  maximum  values.  The  images  are  replotted  in  Figures  52a  to 
52c  after  the  threshold  is  applied  and  we  can  easily  distinguish  each  shape. 
Case  II  demonstrates  similar  results  where  we  applied  a  t i 1  ter  with  lower 
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frequency  cutoff  at  the  15-th  sample  point  and  the  upper  cutoff  at  the 
30-th  sample  point.  The  results  before  and  after  the  threshold  filter  are 
seen  in  Figures  53  to  56c. 


We  have  seen  that  the  radar  restriction  in  frequency  does  not 
eliminate  the  possibility  of  identification.  In  fact  target  images  can  be 
restored  based  on  knowledge  of  the  boundaries. 

Aspect  Windowing 

In  our  application  we  not  only  are  restricted  in  frequency  but  we  are 
restricted  to  measurements  over  a  finite  range  of  aspects.  Therefore  we 
have  applied  a  window  in  aspect  angle,  as  seen  in  Figure  57,  that  is  1 
within  a  specified  angle  and  0  outside  and  is  symmetric  about  the  origin. 
This  rigorously  corresponds  to  data  for  two  aspects  that  are  col  inear  and 
separated  by  180°.  It  has  been  shown,  however,  that  the  transient 
response  at  a  single  aspect  may  be  Fourier  transformed  to  give  approximate 
results  in  both  the  "forward-look"  and  "rear-look"  regions  [24].  Thus,  we 
applied  a  60°  window  in  Case  III  where  the  window  is  ±60°  from  the  nose 
on  aspect  and  symmetric  and  the  resulting  images  before  and  after  the 
threshold  filter  are  seen  in  Figures  58  to  61c.  In  Case  IV  we  applied 
a  30°  window  and  the  corresponding  results  are  in  Figures  62  to  65c. 

For  the  60°  window  the  degradation  is  insignificant  and  for  the  30° 
window  the  degradation  varies  for  each  shape.  As  we  earlier  stated  it  is 
reasonable  to  assume  the  triangular  shape  would  suffer  the  greatest  in 
resolution  because  the  lobes  perpendicular  to  the  sides  were  eliminated  in 
the  30°  windowing  whereas  the  spectra  due  to  the  square  or  circular  shapes 
contained  more  information  and  energy  within  the  window.  This  naturally 
tells  us  that  image  reconstruction  is  dependent  on  the  angle  of  the  window 
axis  as  well  as  the  width  of  the  angular  window.  Even  though  we  lose 
contour  information  at  angles  outside  the  window  we  can  infer  the  identity 
of  the  square  and  circular  targets  from  the  clearly  defined  edges  after 
the  threshold  filter  is  applied.  And  to  some  extent  the  triangular  shape 
is  preserved. 

Frequency  and  aspect  windowing 

A  p-space  window  representing  both  frequency  and  angle  windowing  is 
shown  in  Figure  66.  The  result  of  applying  both  windows  in  cascade  is  seen 
in  Figures  67  to  74c  before  and  after  threshold  filtering.  Case  V, 

Figures  67  to  70c  demonstrates  the  application  of  a  10-30  point  frequency 
filter  and  a  60°  window  and  the  results  are  good  enough  to  distinguish  all 
three  shapes.  Case  VI,  Figures  71  to  74c,  yields  the  results  of  a  10-30 
point  frequency  filter  and  a  30°  aspect  angle  window.  Here  the  square  and 
circular  targets  may  be  defined  and  all  three  may  be  distinguished  but  it 
is  apparent  that  the  most  significant  degradation  results  from  the  appli¬ 
cation  of  a  narrow  aspect  window.  This  leads  to  the  discussion  of  possibly 
extending  the  aspect  angle  window  so  as  to  yield  better  results. 
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Figure  49.  Image  of  triangular  scatterer  after  bandpass  filtering: 
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Figure  51.  Image  of  circular  scatterer  after  bandpass  filtering: 

Flow  =  samP^e  P01’nt  "  ^hiqh  =  sample  point 
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Figure  52c.  Image  of  circular  scatterer  after  bandpass 
filtering  and  threshold  filtering: 

Flow  =  1{)th  samPle  P°int  "  Fhigh  =  3°th  samPle  P°int 
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Figure  53.  Image  of  triangular  scatterer  after  bandpass  filtering: 
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Figure  56c.  Image  of  circular  scatterer  after  bandpass 
filtering  and  threshold  filtering: 

P^,  =  15th  sample  point  -  P".  .  .  =  30th  sai 
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Figure  58.  Image  of  triangular  scatterer  after  aspect  angle  filtering: 
angle  =  60°. 
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Figure  61c.  Image  of  circular  scatterer  after  aspect  angle 
filtering  and  threshold  filtering:  angle=60°. 
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Figure  62.  Image  of  triangular  scatterer  after  aspect  angle  filtering: 
angle  *  30°. 
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Figure  64.  Image  of  circular  scatterer  after  aspect  angle  filtering: 
angle  =  30°. 
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Figure  65c.  Image  of  circular  scatterer  after  aspect  angle 
filtering  and  threshold  filtering:  angle=30°. 
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Figure  67. 


Image  of  triangular  scatterer  after  bandpass  and  aspect  angl 
filtering:  V  =10th  sample  point,  P.  .  .=30th  sample  point, 
angle=60°  ’w  h'gh 
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Figure  68.  Image  of  square  scatterer  after  bandpass  and  aspect  angle 

filtering:  P  =10th  sample  point,  F  =30th  sample  point, 
angle=60° .  low  hl9h 
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Figure  69. 


Image  of  circular  scatterer  after  bandpass  and 
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Figure  70c.  Image  of  circular  scatterer  after  bandpass,  aspect 
angle  and  threshold  filtering: 

^low  *  loth  samPle  P°int  -  ^high  =  3°th  samPle  P°int 
angle  =60°. 
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Figure  71 . 


Image  of  triangular  scatterer  after 
filtering:  F  =10tf1  sample  point, 
angle=30°.  1ow 


bandpass  and  aspect  angle 
^high=3°th  samP^e  P01’nt» 
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Figure  73.  Image  of  circular  scatterer  after  bandpass  and  aspect  angle 
filtering:  P  =10th  sample  point,  F  =  30th  sample  point, 
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Figure  74c.  Image  of  circular  scatterer  after  bandpass,  aspect 
angle  and  threshold  filtering: 

F1qw  =  10th  sample  point  -  Fhigh  =  30th  sample  point 
angle  =  30°. 


F.  Extension  of  the  Aspect  Window 

Since  ease  of  target  identification  is  greatly  facilitated  by 
enlarging  the  viewing  angle  we  would  like  to  find  a  means  for  doing  so 
by  approximate  analysis  techniques.  There  may  exist  several  techniques 
to  arrive  at  a  solution  but  in  this  section  we  have  restricted  ourselves 
to  a  study  of  the  target's  angular  variations. 

Each  particular  scatterer  geometry  has  its  own  characteristic 
spectral  angular  variations  and  because  the  backscatter  spectrum  can  be 
defined  as  the  inverse  Fourier  transform  of  the  target  characteristic 
function  we  can  analytically  calculate  the  backscatter  spectrum  for 
simple  shapes.  The  utility,  of  course,  is  that  we  can  observe  the  angular 
scattering  variations  over  a  limited  range  of  angles  and  by  comparison 
with  the  calculated  variations  we  can  select  the  appropriate  analytic 
function  and  extend  the  window. 


For  the  triangular  target  of  height  a  and  base  width  of  2b  the 
calculated  scattering  function  in  polar  coordinates  is 


r(Pr,e) 


2ab 

T7 
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sin(bP  cose) 

cos(aPrsine)-cos(bPrcose)+iasine^ — tc^ie" 
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"2  2  T 

b“cos  e  -  a  sin  e 


sin(aPrsine) 

aslne 


Jl 


(20) 


and  for  the  square  target  of  height  a  and  width  b  centered  on  the  origin  the 
calculated  scattering  function  in  polar  coordinates  is 


r(Pr,e) 


absi 


a P  cose 
r 


/bP  sine\ 

sin(-ir—  ) 

b^s'ine 
—2 — 


(21) 


and  for  the  circular  target  of  radius  a  centered  on  the  origin  the  calcu¬ 
lated  scattering  function  in  polar  coordinates  is 


r(Pr,e) 


2ira  J-|(aPr) 


(22) 


This  last  one  is  independent  of  e  which  we  expect  since  there  are  no 
angular  variations  in  the  geometry.  Thus  continuation  of  the  angular 
window  for  a  body  of  revolution  would  be  a  matter  of  continuing  the  constant 
value  at  that  frequency.  For  more  complex  structures  we  must  either  calculate 
the  analytic  expression  or  lump  the  results  of  simpler  geometries  to  get 
an  approximate  result.  This  work  is  yet  to  be  completed  but  an  interesting 
sidelight  has  emerged.  In  Figure  75  we  have  plotted  the  angular  variations 
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for  the  triangle  and  square  at  a  fixed  radius  in  P  space.  (The  circular 
shape  would  yield  a  constant.)  The  difference  in  angular  variation  alone 
is  sufficient  to  identify  these  two  targets,  and  to  identify  the  orientation 
of  straight  edges.  Also  the  lengths  of  the  edges  appear  to  be  proportional 
to  the  magnitude  of  the  lobes  and  inversely  proportional  to  their  widths, 
e.g.,  for  an  infinite  planar  surface  the  2-D  scattering  would  be  a  delta 
function. 

F.  Conclusion 

The  studies  and  examples  of  this  section  have  served  to  demonstrate 
some  significant  principles  in  identifying  various  scatterers.  The  first 
principle  is  that  even  though  an  infinite  set  of  frequency  and  aspect  infor¬ 
mation  is  required  in  theory  to  exactly  image  a  scatterer  we  can  produce 
satisfactory  images  with  only  a  finite  portion  of  the  spectrum.  This  is 
true  because  in  the  physical  optics  approximation  the  surface  of  the  scatterer 
gives  rise  to  high  frequency  fields  which  are  preserved  when  a  bandpass 
window  is  applied.  Thus  we  can  determine  the  contour  of  the  object  by 
knowledge  of  the  surface.  The  most  significant  degradation  occurs  in 
constraining  aspect  information  but  we  could  distinguish  the  three  shapes 
even  for  a  30°  window. 

Also,  we  demonstrated  that  the  angular  variations  are  periodic  functions 
in  360°  and  these  yield  insight  into  the  geometry  of  the  scatterer  and 
suggest  the  possibility  of  extending  the  aspect  window  which  could  lead  to 
better  imaging. 

Another  principle  is  that  by  knowledge  of  the  effects  of  windowing 
we  can  determine  the  features  which  are  most  significant,  viz.  the  target 
surface  scattering  etc.  and  thereby  not  try  to  search  for  information  that 
may  be  obscured  by  our  choice  of  window.  We  will  continue  our  efforts  in 
the  areas  of  tapering  the  windows,  extending  the  windows,  perfecting  the 
threshold  filtering,  studying  the  relation  of  geometry  to  the  field  variations, 
filling  in  gaps  in  images  from  knowledge  of  the  known  data  and  determining 
the  limits  of  our  approximations  for  measured  data. 


V.  SUMMARY  AND  CONCLUSIONS 

Research  has  been  done  in  three  major  areas  during  the  first  12  months 
of  this  contract:  Transient  radar  signature  determination  for  cone-like 
targets,  ramp  response  generation  and  limiting  surface  imaging,  and  physical 
optics  inverse  scattering.  None  of  these  efforts  have  been  completed, 
but  some  preliminary  results  are  available  in  each  area. 

Complex  radar  cross-section  data  have  been  obtained  for  six  cone-like 
target  shapes  over  a  10:1  frequency  band,  where  cone  length  £  lies  in  the 
wavelength  range,  A/6  <_  *  <_  1.66A.  The  accuracy  and  scope  of  these  data 
are  considered  to  be  an  advance  in  the  state-of-the-art,  and  a  paper 
summarizing  the  results  will  be  submitted  for  sponsor  approval  in  the  near 
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future.  Cross-section  calculations  have  been  made  for  a  subset  of  the  targets, 
using  Rayleigh  theory  in  the  range  A  ^  A/4  and  geometrical  theory  of 
diffraction  and/or  equivalent  current  approximation  in  the  range  a  >_  A. 

Both  measured  and  calculated  cross-section  amplitude  data  have  been 
compared  to  independently  measured  results  for  cone  shapes  nearly  identical 
to  members  of  our  target  set.  Based  on  all  results  and  comparisons  to  date, 
we  have  more  confidence  in  the  calculated  data  at  a  =  x/6,  where  the 
targets  had  extremely  low  cross-section  for  accurate  measurements,  but  where 
the  Rayleigh  theory  has  been  shown  in  the  past  to  be  accurate.  Above  a  =  x/6, 
the  measured  data  appear  to  be  accurate  within  -  1  dB  and  -  10°.  In 
particular,  the  analytic  technique  used  here  are  not  appropriate  in  the  range 
A/2  <  i  <X  although  the  measured  data  seem  to  smoothly  join  the  calculated 
low-frequency  (Rayleigh)  and  high-frequency  (GTD/equi valent  current) 
predictions. 

The  development  and  calibration  of  the  swept-frequency  system  for 
measurements  in  the  range  1.33X  <.  A  <_  13. 3A  has  been  completed.  Its 
accuracy  is  consistent  with  the  low-frequency  system. 

Time  domain  ramp  response  signatures  have  been  calculated  for  all 
targets  using  the  complex  frequency-domain  cross-section  data.  All  of  the 
waveforms  exhibit  approximate  agreement  with  the  target  profile  function 
(cross-section  vs.  distance)  in  their  forward  portion.  The  latter  part  of 
the  time-domain  waveforms  show  some  indication  of  target  resonance  returns. 
However,  for  the  approximate  nose-on  incidence  angles  of  these  data,  it  is 
obvious  that  cone  resonances  are  not  being  strongly  excited.  The  wave¬ 
forms  demonstrate  that  target  detail  characteristics  (such  as  edge  rounding) 
are  not  resolved,  and  this  is  consistent  considering  the  frequency  range 
of  the  current  data.  However,  the  length  to  base  diameter  ratio  of 
the  different  cones,  and  cone  size  (volume)  can  be  inferred  from  the 
waveforms. 

Some  ramp  response  waveforms  were  also  obtained  using  calculated  data. 
Approximate  agreement  between  experimental  and  calculated  waveforms  is 
observed.  It  appears  that  the  waveforms  from  GTD  calculations  more  nearly 
match  the  measured  data  waveforms  in  the  region  following  the  base  return, 
while  the  base  return  itself  is  more  accurately  approximated  by  the 
equivalent  current  calculation. 

Images  have  beer,  generated  for  the  targets  using  the  limiting  surface 
approach  with  the  ramp  response  signatures  as  input  data.  The  images 
obviously  portray  the  gross  features  of  the  targets.  They  also  show 
surface  ripples  which  are  related  to  the  sharp  cutoff  of  frequency  domain 
data.  A  measure  of  image  accuracy  is  obtained  by  comparing  profile 
function  data  for  the  image  to  the  input  profile  function  data  from  the 
ramp  response  waveforms  at  each  aspect  angle  used  in  the  image  generation 
process.  This  indicator  is  seen  to  correspond  roughly  to  subjective  visual 
evaluation,  although  the  visual  cues  used  in  evaluating  the  image  seem  to 
be  related  far  more  closely  to  silhouette  edge  and  shape  details  rather 
than  any  area  or  volume  parameters. 
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If  the  laser  down-conversion  research  confirms  the  predicted 
characteristics  of  its  electromagnetic  transient  response  return,  then  the 
limiting  surface  imaging  approach  would  be  directly  applicable  to  the 
return  waveforms.  Therefore,  some  research  on  the  basic  limiting  surface 
image  process  has  been  initiated,  concentrating  on  two  areas:  1)  automated 
iterative  image  improvement  algorithms  and  2)  image  confidence  estimates. 

As  the  preliminary  results  in  this  report  show,  the  relationship  between 
image  accuracy  and  the  parameter  values  used  in  image  creation  is 
relatively  well-behaved.  Thus,  automated  iteration  of  these  parameters 
has  been  shown  to  be  feasible. 

It  has  been  shown  in  this  report  that  for  input  data  at  a  single 
look  angle,  several  different  images  can  be  generated  whose  volume  accuracy 
estimates  are  virtually  identical.  This  is  because  the  whole  problem  is 
"underspecified"  in  this  case.  The  concept  of  image  confidence  is  intro¬ 
duced  to  explain  this  situation.  General  analysis  of  image  confidence 
has  been  initiated,  but  has  not  produced  quantitative  results  so  far.  How¬ 
ever,  it  is  concluded  that  qualitatively,  the  image  confidence  is  directly 
related  to  the  number  of  input  waveforms,  and  their  spread  of  frequencies 
and  aspect  angles.  It  is  also  seen  that  the  image  confidence  is  related  to 
the  shape  of  the  target,  in  ways  not  yet  well  defined. 

A  physical  optics  inverse  scattering  theorem  has  been  described.  A 
target  characteristic  function  in  real  three-dimensional  xyz  space  has  been 
shown  to  be  related  to  an  electromagnetic  scattering  characteristic  function 
in  three-dimensional  p-space  (p-space  having  radial  dimensions  from  an 
origin  proportional  to  frequency  and  directional  coordinates  proportional 
to  interrogating  signal  direction  with  respect  to  the  target  in  xyz  space) 
by  a  Fourier  transform.  A  two-dimensional  implementation  of  this  theorem 
has  been  programmed,  along  with  a  computer  intensity  display  algorithm. 
Characteristic  functions  for  some  simple  ideal  two-dimensional  shapes 
have  been  plotted.  Finally,  windowing  functions  have  been  applied  to  the 
p-space  data,  equivalent  to  eliminating  signature  data  in  certain  regions 
of  look  angle  and  frequency.  The  results  demonstrate  which  features  of 
the  target  are  preserved  and  which  are  degraded  as  a  function  of  the 
available  signature  data. 

Advances  to  be  made  in  each  of  these  areas  during  the  final  months 
of  the  contract  should  combine  to  achieve  all  of  the  goals  of  the  contract. 
Approximately  45%  of  the  total  funding  is  available  to  support  the 
remaining  effort. 

Concerning  signature  data,  the  swept  frequency  measurements  shall  be 
completed  giving  data  over  an  80:1  total  bandwidth.  Analytical  results, 
which  become  more  and  more  accurate  for  these  simple  shapes  at  higher 
frequencies,  can  then  be  added  to  the  measured  data  to  obtain  almost  any 
specific  frequency  range  desired.  A  technical  report  summarizing  all  data 
will  probably  be  written. 
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Time  domain  transient  waveforms  will  be  generated  using  the  higher 
frequency  data.  Transient  signature  characteristics  relating  to  target 
resolution  and  target  shape  features  will  be  studied. 

Limiting  surface  images  will  be  generated  using  the  high-frequency 
signature  data.  Also,  research  in  quantitative  image  confidence  esti¬ 
mation,  and  automated  iteration  of  the  process  to  achieve  the  optimum  image 
will  be  concluded.  A  technical  report  and  a  paper  will  be  written  in  this 
area . 

Finally,  the  inverse  scattering  research  will  be  pursued,  buth  to 
gain  insight  into  signature  feature  selection  for  target  identification 
and  techniques  for  optimum  identification  and  imaging  with  limited  signature 
data.  Limited  three-dimensional  images  will  be  produced  using  this  approach 
and  the  best  available  measured  signature  data.  A  technical  report  and  a 
publication  are  expected  to  result  in  this  area. 

The  technical  reports  and  papers  to  be  published,  which  represent 
the  scientific  advances  made  during  this  research,  have  already  been 
discussed.  The  possible  implementation  of  the  results  in  conjunction 
with  a  LASER  radar  system  for  target  imaging  and  identification  has  also 
been  mentioned.  In  addition,  other  areas  of  possible  application  of  the 
work  to  real-world  radar  systems  are  becoming  evident. 

The  data  contained  in  this  report  can  be  applied  directly  to  H.F. 
radar  signature  studies  for  full-size  targets.  Thus,  the  research  is 
applicable  to  target  identification  using  over-the-horizon  radar  systems. 

As  the  results  of  Sections  II  and  III  show,  the  low  frequency  regime 
transient,  signatures  can  be  used  for  target  volume  discrimination  with 
t  10%  accuracy,  and  also  for  determination  of  target  gross  shape  character¬ 
istics.  Other  studies  have  shown  that  the  same  conclusions  apply  to  other 
shapes  such  as  spheres,  cylinders,  cone-cylinders,  etc.  In  fact,  it 
should  be  feasible  to  characterize  the  exhaust  plume  of  an  ascending 
missile  as  well.  In  an  application  where  actual  target  images  are  not 
necessary  but  discrimination  among  known  targets  is  desired,  then  only 
a  small  subset  of  the  return  signature  data  is  needed.  Other  studies 
concerning  aircraft  shapes  have  shown  that  HF  discrimination  is  dependent 
on  the  selection  of  proper  features  for  each  target  class  [25].  The 
imaging  and  inverse  scattering  analysis  of  this  effort  is  directly 
applicable  to  applied  research  on  optimum  feature  selection  for  any  set 
of  targets. 

A  second  area  of  possible  application  makes  use  of  our  experimental 
capability  to  obtain  signature  data  at  arbitrary  bistatic  angles  and 
polarization,  as  well  as  throughout  the  radar  freouency  range.  Future 
measurements  could  be  made  which  are  directly  applicable  to  a  target 
detection  system  using  one  or  a  few  high-power  transmitters  and  several 
dispersed,  inexpensive  receiver  systems.  Our  scattering  signature  vs.  angle 
analysis  results  in  Section  IV  indicates  that  bistatic  return  vs.  angle 
information  in  one  or  two  narrow  frequency  bands  may  in  fact  be  appropriate 
features  for  rapid  target  discrimination  of  cone-like  targets. 
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APPENDIX  A 

TABULATIONS  OF  AMPLITUDE  AND  PHASE  DATA 


The  amplitude  (cm)  and  phase  (deg)  of  each  radar  backscatter 
measurement  is  given  in  this  appendix.  Each  table  that  follows  gives 
the  data  for  a  single  target  along  with  the  data  from  the  check  sphere 
associated  with  the  data  set.  The  data  on  the  target  is  divided  into 
six  groups  of  ten  measurements,  where  the  heading  of  each  group  of  ten 
gives  the  polarization  (V.P.  for  vertical  polarization  and  H.P.  for 
horizontal  polarization)  and  the  aspect  angle  (deg)  (defined  earlier  in 
this  report).  Each  amplitude  and  phase  value  is  preceeded  by  the  harmonic 
number  for  that  value.  The  phase  convention  is  such  that  items  closer  to 
the  antennas  will  give  larger  phase  values. 
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APPENDIX  B 

RAMP  RESPONSE  WAVEFORMS  FOR  TARGET  SET 


The  following  waveforms  are  ramp  response  signatures  for  each  of 
the  targets  used  in  this  study.  Each  graph  contains  ramp  responses  for 
four  look  angles;  <j>=0,  10,  20  and  30°.  The  dashed  curves  are  the  ramp 
responses  as  measured.  The  solid  curves  have  been  altered  by  setting  the 
first  harmonic  phase  terms  to  180°.  Graphs  are  included  for  vertical  and 
horizontal  polarization. 


j 


Figure  B-8.  Target  4  horizontal  polarization  ramp  responses. 
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